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[^^ Abstract. We have analysed near-infrared NTT/SofI observations of the starburst cluster Westerlund 1, which is among the 

_ most massive young clusters in the Milky Way. A comparison of colour-magnitude diagrams with theoretical main-sequence 

^-^ and pre-main sequence evolutionary tracks yields improved extinction and distance estimates of Aks = 1 . 1 3±0.03 mag and d = 
3.55±0.17kpc (DM = 12.75±0. lOmag). The pre-main sequence population is best fit by aPalla & Stabler isochrone for an age 

I I of 3.2 Myr, while the main sequence population is in agreement with a cluster age of 3 to 5 Myr An analysis of the structural 

f^ parameters of the cluster yields that the half-mass radius of the cluster population increases towards lower mass, indicative of 

^L| the presence of mass segregation. The cluster is clearly elongated with an eccentricity of 0.20 for stars with masses between 10 

JL and 32Mq, and 0.15 for stars with masses in the range 3 to IOMq. We derive the slope of the stellar mass function for stars 

^ with masses between 3.4 and 27 Mq. In an annulus with radii between 0.75 and 1.5 pc from the cluster centre, we obtain a slope 

"^ of r = -1.3. Closer in, the mass function of Westerlund 1 is shallower with F = -0.6. The extrapolation of the mass function 

j^ for stars with masses from 0.08 to 120 Mg yields an initial total stellar mass of a52,000 Mq, and a present-day mass of 20,000 

' ' to 45,000 Mq (about 10 times the stellar mass of the Orion Nebula Cluster, and 2 to 4 times the mass of the NGC 3603 young 

, cluster), indicating that Westerlund 1 is the most massive starburst cluster identified to date in the Milky Way. 

> 
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' 1. Introduction At extragalactic distances, however, individual cluster 

■f—i members are often unresolved, and physical properties of SSCs 

t^ Clustered star formation is the dominant mode of star forma- _ such as their initial mass function (IMF) and total mass, de- 

O tion in the Universe with about 80% of all st ars originating gree of binarity, or density profile - must be estimated from 

^ in giant molecular clouds (Lada & Ladai2003). The globular their integrated spectra and photometry alone, which in turn 

*k> cluster systems in the Milky Way or M31 might be remnants are dominated by the most massive, luminous stars present. 

/S of an epoch of Super Star Cluster (SSC) formation early in the The contribution of lower mass stars to the total cluster mass, 

j^ history of these galaxies. Low-mass stars lost in the course of which is crucial for the dynamical evolution and the long-term 

the dynamical evolution of the SSCs could explain the ma- survival of a cluster (e.g., de Grijs & Parmentier 2007 , and ref- 

jority of Pop II field stars found in galaxy halos (Ki'oupa & erences therein), remains uncharted. Indirect estimates of the 



Boily|2002). While the Galactic globular clusters formed more total cluster mass can be derived by extrapolating the IMF to- 

than 12 Gyr ago, 100s of young SSCs with masses compara- wards lower masses. Another possibility is to determine the dy- 

ble to or greater than Galactic globular clusters; i.e., present- namical mass of a cluster from its half mass radius and velocity 

day masses > 100, 000 Mq, have been identified in interact- dispersion under the assumption of virial equilibrium. Studies 

ing galaxi es Uke , e.g., the Antennae galaxies (Whitmore & by different groups either agreed with a Kroupa ( |2002l l IMF 

Schweizer[1995]l. over the whole range of stellar masses (e.g. Larsen & Richtler 

120041) or found a truncated IMF towards lower stellar masses 

(e.g., Rieke et al l 1993 1 McCrady et al. 2005). Young, massive 

Send ojfprint requests to: Wolfgang Brandner clusters, however, are not necessarily in virial equilibrium, as 

* Based on observations collected at the European Southern u . ■ • »• a ■ a t r-. * ■ ti, • .. -ji 

^ photoionization and winds from O-stars in their centres rapidly 
Observatory, La Silla, Chile, and retrieved from the ESO archive (Prog 

ID67.C-0514). 
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expel any remnant gas left over from the formation process 
(Kroupa & Boilv 2002 ; Bastian & Goodwin 2006) . 

The determination of the IMF slope as well as the poten- 
tial presence or absence of a low mass stellar component, is 
critical for constraining the long term evolution of such clus- 
ters. This addresses the question of whether (proto-) Globular 
Clusters are still forming in the Universe today. Moreover, as 
the majority of low mass field stars might have formed in clus- 
ters, the uncertainty in the initial ratio of high to low mass stars 
in clusters results in a corresponding uncertainty in the star for- 
mation and chemical enrichment history of galaxies like our 
Milky Way. 

Therefore we examine starburst clusters in the Milky Way 
and its companion galaxies, where SSCs can be resolved into 
individual stars, and cluster properties like the IMF can be de- 
termined directly. As a starburst cluster we consider any clus- 
ter which is massive enough to house several early O-type stars 
with initial masses in excess of 50 Mq in its centre. Among 
the most massive Galactic starburst clusters studied to date are 
NGC 3603YC (e.g., Brandl et al. T^W Stolte et al. "2006^ and 
Westerlund 2 (Wd 2, Ascenso et al. 2007), both in the Carina 
spiral arm, and the two massive clusters Arches (Stolte et al. 
1^021 120051 Figer et al. l2002l Kim et al. l2006] l and Quintuplet 
(e.g., Cotera et al. 119961 Figer et al. l2004l l in the Galactic centre 
(GC) region. 

One questions of particular interest is whether or not the 
presence of high mass stars hinders or even quenches the for- 
mation of low mass stars. Intense UV radiation fields and 
strong stellar winds of nearby OB stars might lead to the de- 
struction of protostellar cores before low-mass stars can form. 
Furthermore, in extreme environments like the GC region, vari- 
ations in the physical conditions leading to higher thermal 
and/or turbulent pressure, as well as the presence of strong 
magnetic fields might preferentially yield higher mass prestel- 
lar cores than in more benign environs (Morris 119931 Stolte et 
al. 12005, Larson |2006l Klessen et al. 2006) . 

Along with the Arches and Quintuplet clusters in the GC 
region, and the NGC 3603YC in the Sagittarius-Carina spiral 
arm, Westerlund 1 (Wd 1), located in the Scutum-Crux spiral 
arm, is one of the very few young, massive starburst clusters in 
the Milky Way. 

Discovered by Westerlund (1961^ subsequent optical ob- 
servations of the galactic open cluster Westerlund 1 (Wd 1) 
suggested an unusually rich population of cool and hot su- 
pergiants (Westerlund 1987), from which a large cluster mass 
could be inferred. Recent optical spectroscopic and photo- 
metric observations have confirmed these assertions (Clark 
& Negueruela i2002l Clark et al. '200? (henceforth C05); 
Negueruela & Clark l2005T l. Located close to the Galactic plane 
(b = -0.35°) at a distance of ~3-5 kpc, Wd 1 is subject to sig- 
nificant foreground extinction. Thus optical studies of the clus- 
ter were limited to the most massive stars of Wdl. C05 spectro- 
scopically identified about 50 cluster members with masses in 
excess of 30 Mq. Adding 150 probable cluster members identi- 
fied by means of photometry, they estimate that the high-mass 
stellar content of Wd 1 alone amounts to 6000 M©. Assuming a 
Kroupa IMF to extrapolate from this population, C05 then infer 



a likely total mass of ~ 10^ M© for Wdl; directly comparable 
to the masses of SSCs observed in other galaxies. 

Up to now, however, only the evolved, high-mass stellar 
population of the cluster has been characterised. The presence 
of Wolf-Rayet stars (e.g., Crowther et al. 120061 Skinner et al. 
2006a), a pulsar indicative of a recent supernova (Skinner et al. 
2006b, Muno et al. '2006a'), as well as O supergiants, suggest 
an age between 3 and 5 Myr for the cluster. Because of a visual 
foreground extinction of Ay ~ 10 mag, and the uncertainties in 
the intrinsic luminosity of the evolved, massive stars, distance 
determinations to Wd 1 range from 2 to 5.5 kpc (C05). 

Because of its proximity, Wdl provides the unique oppor- 
tunity to directly determine the spatially resolved IMF of a 
SSC and, given that it is located well away from the GC in a 
supposedly 'quiescent' region of the Galactic disc, potentially 
study physical properties imprinted by the formation process 
from the parental molecular cloud on the cluster. With the aim 
to better determine basic cluster properties such as age, fore- 
ground extinction and distance, to directly measure the IMF, 
and to study the structural properties of Wd 1, we have anal- 
ysed deep near IR observations. This paper describes the anal- 
ysis of NTT/SofI observations of Wdl covering a ~ 5' x 5' 
field of view centred on the cluster core, which are sensitive to 
~ solar mass stars at a distance of 4 kpc; a second paper details 
subsequent "VLT/NAOS-CONICA adaptive optics observations 
of selected core regions sensitive to ~0.2Mq. 

The structure of the paper is as follows: in section 2 we 
present the observations and the data analysis. The basic astro- 
physical quantities extinction, distance, and age of Wd 1 are 
derived in section 3. Section 4 discusses the mass function and 
the total stellar mass of Wd 1. Dynamical and structural pa- 
rameters of Wd 1 are addressed in section 5, followed by a 
comparison of Wd 1 to other well studied starburst clusters in 
section 6. Section 7 summarises the results. 

2. Observations and data reduction 

2.1. NTT/SofI observation 

JHKs-broad band imaging observations of Wd 1 (centred on 
RA(2000) = 16''47'"03', DEC(2000) = -45°50'37") and a 
nearby comparison field (located x T to the east, and s; 13' 
to the south of Wd 1, and centred on RA(2000) = 16''47'"43^ 
DEC(2000) = -46°03'47"), each covering an area of 4.5' x 4.5' 
(corresponding to 4.8 pc x 4.8 pc at the distance of Wd 1), 
obtained with NTT/SofI (PI: J. Alves) were retrieved from 
the ESO archive (see Fig. [T]). SofI is equipped with a Hawaii 
HgCdTe detector. The observations were obtained with a plate 
scale of 0.29"/pixel. Individual integration times were 1.2s 
(DIT) and 10 frames were co-added, resulting in 12s of inte- 
gration time per dither position. For each field and filter a total 
of 10 dither positions were used, resulting in a total integration 
time of 120s in each filter. 



2.2. Eclipse raw data reduction 

Data reduction was carried out using the eclipse jitter routines 
(Devillard i2001| l. Because of the high degree of crowding in 
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Fig. 1. Composite colour JHKs images of Westerlund 1 (left) and the comparison field (right) obtained with NTT/Sofl. The field 
of view is 4' x 4'. North is up, and east is to the left. 



the cluster field, sky frames derived from the comparison field 
were used in the reduction of the Wd 1 frames. The resolution 
on the final images is a; 0.75" to 0.80" (see Table [T] for more 
details). 



Table 1. Observing log for the NIR data. The observations were 
caiiied out with the ESO NTT and Sofl on 9 June 2001. 
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seeing 


sat. limit 
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[mag] 


Westerlund 1 


J 


120 


0.81 


10.1 


Westerlund 1 
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120 


0.82 


9.5 


Westerlund 1 


Ks 


120 


0.73 


9.1 


Off-Field 


J 


120 


0.77 


10.2 


Oif-Field 


H 


120 


0.77 


9.6 


Off-Field 


Ks 


120 


0.74 


9.1 



2.3. Photometric analysis and calibration 

PSF fitting photometry in J, H, and Ks for 7000 stars in the Wdl 
area, and 5300 stars in the off-field was derived using the IRAF 
implementation of DAOPHOT (Stetson [T987] l. Photometric ze- 
ropoints and colour terms were computed by comparison of 
instrumental magnitudes of relatively isolated, bright sources 
with photometry from the 2MASS point source catalogue. 

The resulting Ks vs. J-Ks colour-magnitude diagrams 
(CMD) for Wdl and the off-field for all stars with DAOPHOT 
fitting errors <0.2mag are presented in Fig. [2] The dashed line 
indicates the saturation limit (see Table[T]i. The off-field is char- 
acterised by a blue main-sequence of foreground stars, and a 
red sequence of background stars, most of which could be red 



giants located in the galactic bulge. The Wd 1 field shows the 
main-sequence and pre-main sequence members of the cluster 
in addition to the two populations visible in the off-field. Note 
that the red background population appears to be less reddened 
close to Wd 1 compared to the off-field. This is most likely due 
to varying extinction along the line of sight towards the bulge 
at distances larger than the distance to Wd 1 . 

In the following analysis, only stars with DAOPHOT pho- 
tometric fitting eiTors of 0.2 mag or less are considered. 

2.4. Incompleteness characterisation 

Incompleteness simulations with the aim to assess the effect 
of crowding on our abiUty to detect faint sources were car- 
ried out both for the jitter combined off-field and the cluster 
field. In order not to change the crowding characteristics of the 
frames, for each run only 50 stars were added using addstars 
under DAOPHOT. Artificial star magnitudes were allowed to 
scatter ±0.5 mag around a predefined value, and the stars were 
placed at random positions on the Ks-band frames. The detec- 
tion experiment was following the same steps and using the 
same PSF as the initial DAOPHOT analysis. Only those stars, 
whose recovered magnitudes match to within +0.5 mag the in- 
put magnitude, were counted in the final analysis. For each 
predefined magnitude, 10 such frames were created and anal- 
ysed, i.e., 10x50 stars - 500 stars/mag bin. Then the magni- 
tude was increased by 1.0 mag, and 10 new frames were cre- 
ated and analysed. This procedure was repeated 6x, probing 
the magnitude range Ks = 1 1 .9 mag to 17.9 mag for the cluster 
frame, and Ks = 13.9 to 19.9 mag for the field frame. For the J- 
band incompleteness simulations, stars were placed at the same 
positions as on the Ks-band frames. J-band magnitudes were 
computed assuming J-Ks = 2.0 mag, which is an intermediate 
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Fig. 2. Colour magnitude diagram of the cluster (including field contamination, left), and the comparison field (right). The dashed 
line marks the saturation limit. 



colour between the (lower) main-sequence stars in Westerlund 
1 with J-Ks =s 1 .5 mag, and the pre-main sequence stars with 
J-Ks ^2.3 mag. 

In total 240 frames with artificial stars added were anal- 
ysed. The results are summarised in Fig.l3] Note that the lack of 
bright stars in the off-field pushes the 50% completeness limit 
to fainter magnitudes compared to the cluster field. 

2.5. Statistical field subtraction 

With Galactic coordinates 1 * 339.55° and b ^ -0.40°, Wd 1 is 
embedded in a rich population of fore- and background stars. 
The more than 4000 stars detected in J and Ks in the compari- 
son field with DAOPHOT fitting eiTors < 0.2 mag (Fig.|2] right) 
give a good estimate of the field star population mix. This pro- 
vides a firm basis to remove the field star contamination in the 
cluster CMD, and hence determine a clean, relatively unbiased 
cluster population for further analysis. 

The statistical field subtraction is based on a comparison 
of the cluster and the field CMD. The CMDs are subdivided 
into giid cells with a step size of 0.5 mag in colour and mag- 
nitude. The number of field stars within each cell is counted, 
normalised to the ratio of the sky areas covered by the image 
and the selected cluster annuli, and corrected for the differ- 



ences in completeness fraction between the field and the cluster 
(see Fig. [3]l. Finally, the same number of stars is subtracted at 
random from the coiTesponding grid cell in the cluster CMD. 
Examples of the cleaned-up ("field subtracted") cluster CMDs 
are shown in Fig. |4] 

3. Extinction, cluster age, distance and 
evolutionary phases of cluster members 

A starburst cluster most likely forms in a single event ("burst"), 
resulting in a small age spread of its stellar population. Galactic 
starburst clusters like NGC 3603 or Arches with several 1,000 
to 10,000 stars spanning a factor of 100 or more in individ- 
ual stellar mass have the advantage that global properties like 
foreground extinction, distance or age can be derived with 
a relatively small uncertainty by comparison to theoretical 
isochrones. For a detailed analysis of the mass function and the 
dynamical properties of Wd 1, these basic astrophysical quan- 
tities have to be derived first. 

For the following analysis we use Geneva isochrones 
(Lejeune & Schaerer 1200 11 1 describing main-sequence and 
post-main sequence evolution of stars with masses between 0.8 
and 120 Mo, and Palla & Stabler isochrones (Palla & Stabler 
119991 ) describing pre-main sequence and main-sequence evolu- 
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Fig. 3. Results of the completeness simulations for Westerlund 1 . The two figures on the left show the input magnitudes and 
colours and the recovered magnitudes and colours, respectively, for the cluster frame. From this we also derive the photometric 
uncertainties. The figure on the right shows the 50% completeness limits. For the cluster, increasing crowding towards its center 
pushes the completeness limit (solid line) towards brighter magnitudes. The homogeneous stellar density in the comparison field 
results in an overall 50% completeness level down to Ks = 17.55 mag (dashed line). 



tion of stars with masses between 0. 1 and 6 Mq. In the follow- 
ing paper describing the analysis of the adaptive optics data, 
we will provide a more in-depth comparison of evolutionary 
tracks. Isochrones by Palla & Stahler are preferred over other, 
more recent tracks (e.g., Siess et aI. I2000l l as they provide a bet- 
ter fit to the pre-main sequence / main sequence transition re- 
gion in near infrared CMDs of starburst clusters. In their com- 
parison of NIR observations of NGC 3603 YC with PS99 and 
Yonsei-Yale (Yi et al. l2001l ) isochrones, Stolte et al. (2005) also 
found a better fit of the PS99 tracks to the observed pre-main 
sequence / main sequence transition region. Starburst clusters, 
i.e. clusters housing early O-type stars in their centres, might be 
special in this respect, as the intense UV radiation of the most 
massive cluster members rapidly photo-evaporates any remain- 
ing accretion disks or envelopes around lower mass stars in the 
cluster (e.g., Brandner et al. l2000l and references therein). Thus 
pre-main sequence evolution in starburst clusters seems to be 
well described by non-accreting tracks. For a thorough discus- 
sion of short-comings of theoretical isochrones compared to 
observations we refer to Mayne et al. ( I2007I I. Using the same 
set of isochrones as previously used for the analysis of the 



mass functions of NGC 3603 and Arches also facilitates a di- 
rect comparison of the 3 clusters. We use a mass of 6.0 M© as 
the transition point from Geneva to Palla & Stahler isochrones 
as both sets of tracks predict the same NIR magnitudes for stars 
of this mass in the age range considered. 



Table 2. Basic astrophysical parameters extinction, DM (dis- 
tance), metallicity and age of Wd 1 as derived from comparison 
of the Ks vs. J-Ks CMD with theoretical Geneva and Palla & 
Stahler isochrones. 



Ak. 
mag 



DM 

[mag] 



age 

Myr 



1.13±0.03 



12.75±0.10 
(3.55±0.17kpc) 



0.02 3.6±0.7 



The values for extinction, distance modulus and age of the 
cluster, which we present in the following, were derived in an 
iterative process by fitting isochrones to the Sofl and NACO 
near-infrared data. 



Wolfgang Brandner et al.: Intermediate to low-mass stellar content of Westerlund 1 



10 k 



12 



W 14 



16 



18 



Wdl after field sub: 1751 stars 

0.0" < rad < 84.7" (O.OOpc < rad < 1.50pcJ 



- r: 

■■■6.' ■ 




Geneva 3.9 Myr 
■ PS99 3,2Myr 



DM = 12.75 mag 
A_Ks = ] . 1 3 inag 



2 3 

J-Ks [mag] 



4.0 Mo 
3.SMa 
2.0 Mo 



1.5 Ma 
1.2 Mo 



Wdl after field sub: 1299 stars 

84.7" < rad < 119.8" (l.SOpc <rad < 2.12pc) 



10 k 



12- 



^ 14- 



16- 




18- 



Geneva '3.^ Myr" 
PSQg 3.2 Myr 



DM = 12.75 mag 
A_Ks = ].13 mag 



2 3 

J-Ks [mag] 



4.0 Mo 
3.3 Mo 
2.0 Mo 



1.5 Mo 
1.2 Mo 



. Wdl after field sub: 760 stars 
119.8"<rad< 188.5" (2.12pe<rad< 3.34pc) 



10 k 



12- 



til 14- 



16- 




18- 



Geiieva.3.9 Myr 
PS99 3.2Myr ' 



DM = 12.75 mag 
A_Ks = 1.13 mag 



2 3 

J-Ks [mag] 



4.0 Mo 

. 3.SM0 

2.0 Mo 



1.5 Mti 
1.2 Mo 



Fig. 4. Field subtracted CMDs of the cluster for 3 annuli. A 3.9 Myr Geneva (Z=0.020, with overshooting and OPAL opacities) 
and a 3.2 Myr Palla & Stahler isochrone are overplotted for comparison. The dashed line marks the saturation limit. Photometric 
uncertainties in K-magnitude and J-K colour as determined by the artificial star tests to determine incompleteness are plotted to 
the left of the cluster sequence. The dotted lines indicate stellar masses of cluster members according to the Geneva isochrone 
for masses > 6.OM0, and the Palla & Stahler isochrone for masses less than 6.OM0. 



3.1. Extinction 

As the intrinsic J-Ks colours of main sequence stars with 
masses in the range 6 to 30 M© just vary between -0.1 and 
-0.2 mag, the foreground extinction can be derived by sim- 
ply fitting a zero-age main-sequence (ZAMS) to the main- 
sequence population of Wdland assuming a standard Rieke 
& Lebofsky (1985) extinction law. In Fig. |5]only stars in the 
Sofl data set with DAOPHOT fitting errors less equal 0.05 mag 
are shown. The artificial star tests carried out as part of the 
incompleteness simulation indicate that the DAOPHOT fitting 
errors underestimate the true photometric errors. Interestingly, 
the recovered magnitudes are on average brighter than the in- 
put magnitudes of the artificial stars. This asymmetry in the 
photometric errors around zero is also present in the colour es- 
timate J-Ks, though less pronounced. The photometric errors 
explain part of the observed scatter in the colour of the main 
sequence stars. By comparing CMDs for different regions in 
our field of view, we also see evidence for differential extinc- 
tion. In general, the regions to the west and south of the cluster 
centre suffer slightly lower foreground extinction than the re- 
gions to the east and north of the cluster centre. As a best fit, 
we get Aks - 1.13+0.03 mag. Since all our measurements are 
in the near-infrared, and the theoretical isochrones used have 
been transformed to JHK-magnitudes and colours, the results 
are relatively insensitive to any deviation from a standard ex- 
tinction law. 

In their analysis of the NIR colours of 1 8 W-R stars in Wd 1 
Crowther et al. (I2006I I derived Aks - 0.96 ± 0.14 mag, which 



within the quoted uncertainties overlaps with the Ks-band ex- 
tinction derived by us. We note, however, that for the two W-R 
stars in overlap with our sample (see Table |4]), Crowther et al. 
(2006 ) observe J-Ks colours 0.35 to 0.46 mag bluer than we do. 
Variability as reported by Bonanos ("2007 ) might explain part of 
this discrepancy, though there might also be a systematic offset 
in the zeropoint calibrations of the two data sets. 

For comparison to other extinction determinations, which 
were exclusively based on studies in the visual, we also com- 
pute the visual extinction. Assuming a normal extinction law 
with R = Ay / Eb-v = 3.1, Aks - 1.13 mag corresponds to 
Av * 10.1 mag (Rieke & Lebofsky fT985T l. 

Westerlund ( |1961| l derived an optical colour V-I - 
+4.5 mag for the ^^0 brightest stars of the cluster. By assum- 
ing an intrinsic colour of (V-I)o = -0.35, which is typical for 
early type stars, this corresponds to E(V-I) - 4.85 mag. If one 
further assumes a standard extinction law with Ay = 1.93 E(V- 
I) (see Rieke & Lebofsky 1985), this results in Ay « 9.4 mag. 
Note that Westerlund ( 119611 ) deduced a higher Ay of 1 1 .2 - 
12 mag by assuming a ratio between total and selective absorp- 
tion of 2. 3± 1.4 (Kron & Mavall [79601 1. Follow-up photometric 
and spectroscopic observations considering about 260 cluster 
members yielded Ay s; 9.7+0.8 mag (Westerlund 1987'). Clark 
et al. (2005) estimate Ay - 11.6 to 13.6 mag, and present ev- 
idence for anomalous extinction towards «20 OB supergiants 
in Wd 1. While for normal interstellar extinction the ratio of 
total to selective extinction R = 3.1, R values in the range 1.7 
to 5.5 have been observed in high-mass star forming regions 
(Patriarchi et al. 120011 but see also Hadfield & Crowther |2006l 
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Fig. 5. Colour magnitude diagram of the main sequence popu- 
lation of the cluster for stars with DAOPHOT fitting uncertain- 
ties smaller than 0.05 mag in J and Ks. A Geneva isochrone 
(soUd line, age of 3.9 Myr, DM - 12.75 mag) for Ak = 
1.13 mag is overplotted. Photometric uncertainties (root mean 
square) as determined by the artificial star test are indicated to 
the left of the cluster main sequence. Even though some of the 
stars with colours 0.3 to 0.5 mag redder than the main-sequence 
have been identified as W-R stars, the relatively wide colour 
range is in part due to differential foreground reddening. 



and Lamzin l2006] l. For R = 3.7, Ak = 1.13 ± 0.03 mag coiTe- 
sponds to Ay ~ 12.0 mag, i.e. in the range of the visual extinc- 
tion values determined by Clark et al. ( 120051 ) towards some of 
the evolved stars in Wd 1 . 

3.2. Age 

For the age determination, we take advantage of the fact that 
the width and the shape of the pre-main sequence - main- 
sequence transition region in J-Ks varies with age. This allows 
us to derive the age of the cluster by comparison with theoret- 
ical isochrones. The best fitting Palla & Stahler isochrone are 
for ages of 3.2 to 3.4 Myr Isochrones for an age of 3.0 Myr 
or 3.6 Myr, respectively do not fit as well, in particular when 
taking the adaptive optics data into account. A more detailed 
comparison of different sets of theoretical isochrones will be 
presented in the second paper, which features a detailed analy- 
sis of the cluster population in the range of 0.4 to 3 Mg. We es- 



timate the age of the low-mass stellar population of Wd 1 since 
crossing the "birthline" (see Stahler ITW?]) to 3.3±0.2 Myr. 

At young ages of < 10 Myr, pre-main sequence evolutionary 
tracks originating at the stellar birthline, and main-sequence / 
post-main-sequence evolutionary tracks originating at the zero- 
age-main-sequence are not necessarily on the same absolute 
timescale, and might show systematic offsets. This is exempli- 
fied by the Orion Nebula Cluster (ONC), for which Meynet et 
al. (I1993I I fit the youngest Geneva isochrone in their grid for 
an age of 4 Myr to the upper main sequence of ONC, while the 
canonical age for the low mass stars is <2Myr (Hillenbrand 
fT997l Palla & Stahler [T999T l. 

In the mass and wavelength (colour) range accessible with 
the Sofl near-infrared data, Geneva isochrones (Lejeune & 
Schaerer l2001l l are rather insensitive to age, in particular as 
near-infrared colours like J-Ks always sample the Rayleigh- 
Jeans tail in the spectra of the hot, massive stars. Consequently, 
we cannot use the Geneva isochrones to derive further con- 
straints on the age of Wdl. Based on the observations of 
evolved stars in Wd 1, both Clark et al. (.200 5 ) and Crowther et 
al. ( 120061 ) argue for an age of 4 to 5 Myr In the following anal- 
ysis, we use Geneva isochrones in the age range 3.0 to 5.0 Myr. 
The isochrones are based on the "classical" evolutionary tracks 
including overshoot, and have been computed for a metallicity 
ofZ=0.02. 

As an average age for the low (3.3 ± 0.2 Myr) and high- 
mass (3.9+ 1.0 Myr) stellar population of Wd 1, we adopt truster 
= 3.6 + 0.7 Myr. This also results in consistent values for the 
brightness and colours predicted by both sets of tracks for the 
6M0 stars. 

3.3. Distance 

Once extinction and age have been established, the distance 
to the cluster follows straight from the observed near-infrared 
magnitude of the transition region between pre-main sequence 
and main-sequence. For the 3.2 Myr PS99 isochrone, we get 
DM = 12.75±0.10mag or d = 3.55 ± 0.17 kpc. 

This is in agreement with the study by Westerlund ( 119871 ). 
who derived a distance modulus of DM = 13.6+0.7 mag, as 
well as Clark et al. (2005), who argue for a distance between 
2 and 5.5 kpc, and Crowther et al. (120061 ). who derive DM - 
13.50+0.66 mag based on the analysis of the NIR colours of 18 
W-R stars in Wd 1. More recently, Kothes & Dougherty (20071) 
derived a distance of 3.9±0.7 kpc based on the radial velocity of 
an HI feature likely associated with Wd 1, assuming a standard 
Galactic rotation curve and a distance of 7.6 kpc (Eisenhauer et 
al. 2005, Reid 1993) between the Sun and the GC. 

The results on extinction, age and distance are summarised 
in Table |2l 

3.4. Evolutionary phase of cluster members 

Wd 1 provides a unique, and quite interesting snapshot of stel- 
lar evolution as a function of mass. At an age of a;3.6Myr, 
the most massive cluster members have already left the main- 
sequence. Only stars with Ks > 11 mag, corresponding to 
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an MK-type of B0.5V (m k 15 Mq), are still close to the 
ZAMS, whereas more massive stars already started to evolve 
off the ZAMS. The cluster main-sequence extends down to Ks 
^ 15 mag, corresponding to an MK-type of AOV (m = 3 M©). 
Stars with masses between 2 and 3 Mq are located in the pre- 
main sequence - main-sequence transition region, and stars 
with masses less than 2 M© are still in the pre-main sequence 
phase. 

Massive young clusters like Wd 1 offer the rare opportunity 
to bring evolutionary tracks for higher mass stars, starting with 
the zero-age main sequence, and pre-main sequence evolution- 
ary tracks starting, e.g., with the stellar birthline (Stabler [l983l l 
on the same evolutionary time scale. This is a crucial empiri- 
cal calibration point in developing a consistent theory of stellar 
evolution encompassing all stars from 0.08 M© to I2OM0. 

4. Mass function and cluster mass 

We compute the centre of Wd 1 from the field-subtracted spa- 
tial distribution of cluster members. We find the cluster cen- 
tre to be located at RA(2000) = 16''47'"04.0% DEC(2000) = 
-45°5r04.9". This location is used as a reference for the fol- 
lowing discussion. 

Table 3. Mass function slope F for stars with masses between 
3.4 and 27 Mq. 



Radial Dist. 


3.2 Myr 


3.9 Myr 


4.9 Myr 


3.9 Myr 




z=0.02 


z=0.02 


z=0.02 


z=0.04 


0-0.75pc 


-0.7 


-0.6 


-0.5 


-0.6 


0.75- 1.5 pc 


-1.4 


-1.3 


-1.3 


-1.4 


1.5- 2.1 pc 


-1.7 


-1.6 


-1.5 


-1.6 


2.1- 3.3 pc 


-1.7 


-1.7 


-1.9 


-2.0 



r values are quoted for different annuli, and three Geneva isochrones 
for z=0.020 and for ages of 3.2, 3.9, and 4.9 Myr, respectively. The 
column to the right lists T based on a Geneva isochrone with z=0.040 
and for an age of 3.9 Myr. There is a slight dependence of T on the age 
and the metallicity of the adopted isochrone. 
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Fig. 7. Cumulative mass function for 4 annuli for stars with 
masses between 3.4 and 32 Mq. 



4.1. Mass function and number of cluster members 

Individual masses for cluster members are derived by com- 
parison of theu" location in the CMD with a 3.9 Myr Geneva 
isochrone for stars more massive than 6.0 Mq, and a 3.2 Myr 
PS99 isochrone for stars with masses of 6.0 M© and less. We 
note that in the evolutionary time scale defined by the tracks 
by Palla & Stabler, a star with a mass of 3.5 Mq takes about 
1.4 Myr to cross from the birthline to the zero-age main se- 
quence (ZAMS), while the location of the birthline for a star 
with 6.0 Mq coincides with the ZAMS. No attempt was made 
to correct for intrinsic infrared excess and/or extinction to- 
wards individual stars. This is justified by the fact that a J- 
H vs. H-Ks colour-colour diagram of the cluster members 
indicates that only a small fraction (less than 10%) have an 
intrinsic infrared excess. Stellar masses were assigned based 
on the JHKs-brightness and -colours. For stars above 3.5 Mq, 
only those with J-Ks colours within -0.3 and 4-0.4 mag of 
the isochrone were considered for the mass function. As the 
pre-main sequence - main-sequence transition region marks an 
area with a broader spread in colour, for stars below 3.5 Mq, the 
J-Ks colour range was extended to allow a variation of -0.3 to 
+0.6 mag around the colour indicated by the PS99 isochrone. 

At an age of 3.9 Myr, a star with an initial mass of 30 Mq 
has lost about 3% or 6.5% of its initial mass for a metallicity 
of Z = 0.02 or 0.04, respectively, according to the Geneva evo- 
lutionary models. While this has a minute effect on the mass 
function, we assign initial rather than present-day masses to 
each of the stars. The mass function obtained this way is then 
corrected for incompleteness, and the slope of the mass func- 
tion is computed for the range 3.4 to 27 Mq. 

Mass functions are traditionally visualised as histograms. 
Deriving the mass function slope from a histogram can be prob- 
lematic, as the value of the fitted slope in general varies with 
bin size and locations (Maiz-Apellaniz & Ubeda l2005l l. Slopes 
derived from cumulative mass functions are much more robust 
(e.g., Stolte et al. 12006b . In the following we discuss both stan- 
dard histogram mass functions and cumulative mass functions. 

Figure l6] shows the mass function of the cluster assuming 
solar metallicities for four annuli. The mass function slope F 
has been fitted for the mass range 3.4 to 27 Mq. With F - -0.6, 
the mass function is shallow in the inner aiO.75 pc with respect 
to a standard Salpeter mass function with F = -1.35, getting 
steeper for radial distances from the cluster centre between 0.75 
and 1.5 pc with F = -1.3, F = -1.6 for radii between 1.5 pc 
and 2.1 pc, and F - -1.7 for radii larger than 2.1 pc. Because 
of the small number of high mass stars, as well as uncertainties 
in the statistical subtraction of field stars, the slopes in the two 
outermost annuli have uncertainties of 5% to 10%, whereas the 
slopes in the inner annuli have uncertainties smaller than 5%. 

Table l3] summarises the derived mass function slopes for 
masses between 3.4 and 27 Mq. The lower limit of 3.4 Mq en- 
sures that also in the innermost annulus only stars with masses 
above the 50% completeness limit are included. The upper limit 
has been selected such that stars with 27 Mq are still below our 
photometric saturation limit even for an age of 4.9 Myr. 

The cumulative mass function for the same four annuli is 
shown in Fig. U\ It confirms that the slope of the mass function 
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Fig. 6. Mass functions (solid line) of the cluster for 4 annuli. The dashed line gives the incompleteness corrected mass function. 
The 50% completeness limit is indicated by a vertical dotted line. The mass determination is based on a 3.9 Myr Geneva isochrone 
with Z=0.020 for stars with masses >6.0Mq, and on a 3.2 Myr Palla & Stahler isochrone with Z=0.019 for stars with mass 
<6.0Mq. The slope of the mass function (dash-dotted line) is computed from stars with masses between 3.4 and 27 Mq. 



is getting steeper with increasing distance from the cluster cen- 
tre. The mass function in the annulus with distances between 
0.75 and 1.5 pc, i.e., the annulus centred on the half-mass ra- 
dius (see sectionlSll provides a good fit of the overall mass func- 
tion of Wd 1. The MF slopes in the two outermost annuli (1.5 
to 2.1 pc and 2.1 to 3.4 pc, respectively) is less well defined. 
As discussed already in the case of the histograms, this is in 
part due to the relatively small number of massive stars in the 
outermost annulus, and in part due to uncertainties in the sub- 
traction of field stars towards lower masses. In the following 
we thus assume an average slope of F - -1.6 for the two out- 
ermost annuli. 

Assuming a single power-law initial mass function for the 
cluster in the mass range 0.5 to I2OM0 for each of the an- 
nuli, we can now estimate the total (initial) number of clus- 
ter stars. After correction for incompleteness, 1,385 stars with 
masses between 3.4 and 30 Mq are present in the cluster Below 
0.5 Mq we assume a Kroupa IMF slope with F - -0.3. Tablel6] 
lists the total number of stars in the cluster Including the ^65 
supernovae, which might already have exploded in Wd 1 (see, 
e.g., Muno et al. 2006b , and the estimate below), the extrapola- 
tion indicates that we can expect a total of a: 100,000 stars with 
initial masses between 0.08 to 120 Mg in Wd 1. 



4.2. Cluster mass 

Table HI lists bright stars close to the saturation limit of Ks = 
9. 1 mag in our sample, which also have spectral types deter- 
mined by Clark et al. ( I2005I ). Mass estimates derived by com- 
parison with Geneva isochrones show that these stars indeed 
have masses around 30 Mq, as also determined by Clark et al. 
for the faint end of their sample. This overlap between the two 
samples enables us to study the cluster mass over the entire 
range from our 50% completeness limit up to the highest mass 
stars identified by Clark et al. ( 120051) . 

We determine the total stellar mass of the cluster in two 
slightly different ways. First, we add the mass of each star de- 
tected in the SofI data, correcting this number by the complete- 
ness correction. The completeness correction was done on a 
star-by-star basis for each star detected, taking both the bright- 
ness and the radial distance from the cluster centre into account 
(i.e. for a 6 Mq star located in the very cluster centre, the correc- 
tion factor applied is higher than for a 6 Mq star located near the 
outer edge of the innermost annulus). Table l5] summarises the 
stellar mass in the cluster for three different annuli and stellar 
mass bins, only counting stars with masses between the respec- 
tive 50% completeness limits and 30 Mq. Main sequence stars 
with masses between 3.5 Mq and 30 Mq amount to ail 1300 Mq. 
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Table 4. Stars in the overlap between our sample based on NIR photometry and the optical spectroscopic studies by Clark et al. 
(I2005l l and Crowther et al. ( 2006) . 



W87 


RA(2000) 


DEC(2000) 


Ks 


J-Ks 


M3.2Myr 


Mj.QMyr 


M4.9Myr 


Spt 


# 






[mag] 


[mag] 


[Mo] 


[Mo] 


[Mo] 




60 


16 47 4.12 


-45 5152.3 


9.13±0.01 


1.59+0.01 


36.7 


32.0 


28.0 


09.5Ia - BO.Sla" 


29 


16 47 4.41 


-45 50 40.0 


9.17±0.03 


1.58±0.03 


36.4 


31.7 


27.8 


09.51a - B0.5Ia 




16 47 6.54 


-45 50 39.2 


9.25±0.05 


2.03±0.05 








WR77s: WN6o 


15 


16 47 6.63 


-45 50 29.8 


9.66±0.01 


1.56±0.01 


30.6 


27.6 


24.8 


OB binary/blend 




16 47 7.65 


-45 52 36.0 


9.33±0.02 


1.90±0.02 








WR77sb(0): WN60'' 



identified as a variable star (Bonanos l2007b 



The 1st column lists the identification according to Westerlund ( !1987t . Mass estimates (columns 7-9) are based on Geneva isochrones for ages 
of 3.2, 3.9 and 4.9 Myr, respectively. Spectral types from Clark et al. (OB-stars) and Crowther et al. (W-R stars) are listed in column 10. 

Table 5. Integrated stellar mass for different annuli and stellar mass bins, and resulting total cluster mass including all stars 
detected on the SOFI frame with individual mass estimates up to 30 Mq. 



Radial Dist. 


0.8 -1.5 Mo 


1.5 -3.5 Mo 


3.5-30 Mo 


M,„, 


pc 


[Mo] 


[Mo] 


[Mo] 


[Mo] 


0-0.75 


7 


>210 


4135 


>4345 


0.75-1.5 


>220 


1090 


4930 


>6240 


1.5-2.1 


455 


780 


1825 


3060 


2.1-3.3 


260 


340 


540 


1140 


total mass 


>935 


>2420 


11430 


>14785 


total mass" 








>20800 


(0.8-120 Mo) 











^ including the estimate of a;6000 Mq in stellar mass for stars more massive than 30 Mq (Clark et al. 12005b . 

Mass bins below the 50% completeness limit were not corrected for incompleteness. As a lower limit, only the total mass of stars actually 
detected in these bins is quoted, preceded by a ">". As the field of view is limited to an area of 4.3 pc x 4.3 pc, approximately centred on the 
cluster, only stars located within 3.3 pc of the cluster centre are considered. 



Including all low-mass stars with masses down to 0.8 Mq, the 
total stellar mass in this mass range is at least a; 14700 M©. If we 
add to this the estimate from Clark et al. (2005 ) that stars more 
massive than 30 Mq contribute another 6000 IVIo, the combined 
stellar mass of Westerlund 1 amounts to at least 2 x W^Mq 
(right column). Higher resolution adaptive optics data, which 
we are going to discuss in a forthcoming paper, reveal that the 
cluster mass function extends well below 0.8 M©, adding even 
more stars and mass to the cluster. 

The second way to estimate the total stellar mass is based 
on the mass function slopes, and the observed number of stars 
with masses between 3.4 and 30 M©, taking the completeness 
correction into account. By assuming a single power law mass 
function in each of the cluster annuli for stellar masses between 
0.5 Mq (the peak of the Kroupa IMF) and 120 Mq, and assum- 
ing r - -0.3 for masses between 0.08 and 0.5 Mq, we can 
extrapolate the total number of stars and the total stellar mass. 
For the extrapolations we carried out Monte Carlo simulations 
assuming a randomly populated IMF. The simulations yield a 
total number of stars of x 100,000 (Table [6]l, and a total ini- 
tial stellar mass of k 52, OOOMq (Table |7]) for the cluster, with 
an uncertainty of about ±15%. The Monte Carlo simulations 
also indicate that an uncertainty in the mass function slope in 
the outermost annulus by ±0. 1 results in a change of the total 
cluster mass by ±5%. An additional uncertainty derives from 



the unknown binary properties of the cluster members. Both 
multiplicity fraction and distribution of mass ratios need to be 
known in order to correct any derived mass function and the to- 
tal cluster mass for binary properties (see, e.g., the discussion 
in Stolte et al. ' 20051 l2006t . Our seeing limited Sofi data with a 
resolution of 0.8" (i.e. ^2800 A.U. at the distance of Wd 1) do 
not provide strong constraints on the binary properties, so that 
we cuiTently cannot quantify the effect of unresolved binaries. 

For comparison, C05 estimated a total mass of 100,000 Mq 
by assuming that Wd 1 currently houses 140 stars with ini- 
tial masses >30 Mq, and then extrapolating down to stars with 
masses 0.08 Mq based on a two-part Kroupa (120021) IMF with 
(in our notation) F = -1.3 for stars with masses >0.5 Mq, and 
F = -0.3 for stars less massive. Mengel & Tacconi-Garman 
(I2007I I derive a dynamical mass of 6.3^3'^ x 10"* Mq from the 
velocity dispersion of four red supergiants in Wd 1 . Our esti- 
mate of a total initial cluster mass of 52,000 Mq within a;2.5 pc 
of the cluster centre is in agreement with these studies. As there 
seem to be cluster members beyond 2.5 pc from the cluster cen- 
tre, the actual cluster mass might be even higher This confirms 
that Wd 1 is indeed one of the most massive starburst clusters 
in the Milky Way, and one of the best local SSC analogues. 

The extrapolation based on the observed slope and the sin- 
gle power law mass function also indicates that the cluster 
initially had a!65 stars with masses between 50 and 120Mq. 
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Table 6. Total number of stars in the cluster based on extrapolation of the mass function slope, and scaling according to the 
observed number of stars with masses between 3.4 and 30 M© for three annuli. 



Radial Dist. 
pc 


Tobs 


Nob. 
(3 .4-30 Mo) 


-'-^ compl 

(3.4-30 Mo) 


No.08-0.5 

(0.08-0.5 Mo) 


No.5-30 

(0.5-30 Mo) 


N30-50 
(30-50 Mo) 


N50-120 
(50-120 Mo) 


N,o, 
(0.08-120 Mo) 


0-0.75 
0.75-1.5 

1.5-3.5 


-0.6 
-1.3 
-1.6 


285 
512 

322 


404 
634 
347 


14800 
32700 
36700 


1585 
8180 
7770 


37 
18 
6 


42 
15 
6 


15460 
40910 
44480 


sum 




1119 


1385 


84200 


17530 


61 


63 


101850 



For the estimate in each annulus, a power-law IMF with the observed F is assumed over the mass range from 0.5 to 120 Mo, and F = -0.3 for 
masses between 0.08 and 0.5 Mo. 

Table 7. Current and total initial stellar mass of the cluster based on extrapolation of the mass function slope, and scaling 
according to the observed number of stars with masses between 3.4 and 30 M© for three annuli. 



Radial Dist. 
pc 


Fobs 


Mo.08-0.5Mo 

[Mo] 


Mo.5-30A(o 

[Mo] 


Mi30_50Mo 

[Mo] 


M50-120M0 
[Mo] 


M,o, 
[Mo] 


0-0.75 

0.75-1.5 

1.5-3.5 


-0.6 
-1.3 
-1.6 


3110 
6940 
7750 


5460 
12550 
9400 


1430 
690 

275 


3210 
1080 
430 


13210 
21260 
17550 


sum 




17800 


27410 


2395 


4720 


52320 



For the estimate in each annulus, a power-law IMF with the observed F is assumed over the mass range from 0.5 to 120 Mo, and F = -0.3 for 
masses between 0.08 and 0.5 Mo. 



For a cluster age of x!4Myr, these stars should already have 
turned into supernova. As discussed by Muno et al. (2006b), 
this corresponds to a supernova rate of one every ^ 15,000 yr. 
The main effect of the supernovae on the cluster would have 
been the removal of any remnant gas tracing back to the for- 
mation of the cluster as well as a loss of si4,700 M© in stel- 
lar mass (corresponding to a!lO% of the initial stellar mass in 
the cluster). According to the IMF extrapolations, stars with 
masses between 30 and 50 Mq should have initially contributed 
~2,4OOM0 to the cluster mass. These stars, many of which are 
by now W-R stars, should have lost a total of 1 ,000 to 1 ,500 Mq . 
Thus at the high-mass end, the cluster should have lost a total 
of 5,700 to 6,200 Mq in stellar mass since its formation, result- 
ing in a current cluster mass of ai45,OOOM0. Nevertheless, ad- 
ditional thus far undetected cluster members might be located 
outside the field of view of the Sofl data. Thus the actual mass 
of Wd 1 could be even higher. 

4.3. Effect of age uncertainties and metallicity on 
mass function slope and cluster mass 

As discussed before, the uncertainty in the age determination 
of Wd 1 traces back to the fact that we cannot a priori assume 
the Geneva and Palla & Stabler isochrones to be on the same 
absolute timescale. Based on the width and shape of the pre- 
main sequence / main-sequence transition region, the best fit- 
ting PS99 isochrone can be determined with a small uncertainty 
in age. For stars with masses above 6Mq, the selection of the 
best fitting Geneva isochrone is not as straight forward. 

Crowther et al. (120061 1 discuss that Wd 1 might have higher 
than solar metallicity as it is located at smaller galactocentric 
distances than the Sun and formed just recently. Based on the 
analysis of the integrated spectrum of Wd 1, however, Piatti et 
al. ( I1998D conclude that it has close to solar metallicity. In order 



to get an estimate on the effect of metallicity on the derived 
mass function slope, we also applied a Geneva isochrone for an 
age of 3.9 Myr and twice solar metallicity to the observations. 

Table |4] exemplifies that the selection of different Geneva 
isochrones with ages between ^3 and 5 Myr and solar metal- 
licities as well as the selection of a twice solar metallicity 
isochrone for an age of 3.9 Myr results in a variation of the 
mass estimate by ±15% for the most massive stars in our 
sample. Towards lower masses this effect becomes less pro- 
nounced, and only the upper most bins of the mass function 
are affected. The estimate of the integrated cluster mass is even 
less subject to the uncertainties in the age estimate. Table |3] 
summarises the different mass function slopes F derived by ap- 
plying Geneva isochrones with ages of 3.2, 3.9, and 4.9 Myr, 
respectively, to the data. 

From this we can conclude that the determination of the 
slope of the mass function as well as the cluster mass are quite 
robust, and relatively insensitive to the remaining uncertainties 
in the age and metalUcity of Wd 1 . 

5. Structural and dynamical properties of Wd 1 

5.1. Cluster size and shape 

The identification of several 1000 cluster members with masses 
between 0.8 and 32 Mq facilitates a study of the structural and 
dynamical properties of the cluster. 

From the distance of 3.55 kpc between the Sun and Wd 1, 
a distance of 8.0 kpc between the Sun and the GC, and the 
Galactic coordinates of Wd 1, we compute its galactocentric 
distance to 4.8 kpc. As already pointed out in the introduction, 
Wd 1 is thus at a large enough distance from the GC that tidal 
effects leading to stripping of stars can be neglected. 

Only counting stars with masses between 3.4 and 32 Mq, 
the stellar surface density in the cluster core is at least 3.6 x 
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lO-'Mo/pc^. The surface density profile follows a oc (l + 
{rjaf^y^ law with a - 1.1 pc (see Fig. ^ left). According to 
Elson et al. (I1987I eq. 22) this corresponds to a core radius of 
0.4 pc. We do not attempt to fit a King profile, as the crowd- 
ing and saturation of the brightest stars near the cluster centre 
prevents any accurate determination of the cluster's core radius 
based on the Sofl data. 

The half mass radius for all stars with masses in the range 
3.4 to 32 Mq, averaged over all position angles, is 1 .0 pc (Fig. IS] 
right). Figurel9]gives evidence that the cluster is not spherical. 
Here we compute the half mass radius as a function of position 
angle, taking only stars located either in a quadrant ±45° rela- 
tive to the position angle (PA) or located in the opposite quad- 
rant (PA ±180°) into account. The mass range of 3.0 to 32 Mq 



has been divided into two equal size bins in log mass. In both 
mass bins, the cluster is clearly elongated with its major axis 
along PA a; -10° to 10° and its minor axis along PAx; 80° to 
90°. The elongation is more pronounced for stars with masses 
in the range 10 to 32 Mq with an eccentricity of the cluster e - 
1 - SMoraxis ^ 0.19 (Fig.lgl left). The average half-mass-radius 

majoraxis ^ ^ rj ^ '^ 

for stars in this mass bin is 0.8 pc, which is in perfect agree- 
ment with the half-light radius reported by Mengel & Tacconi- 
Garman (I2007I I based on another set of NTT/SofI observations 
of Wd 1. For stars with masses in the range 3.5 to 10 Mq the 
eccentiicity of the cluster is e = 0.15 and the average half mass 
radius is 1 . 1 pc (Fig. l9] right). The larger half mass radius in the 
lower mass bin indicates the presence of mass segregation. 
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An elongation was also reported by Muno et al. (2006b') in 
their analysis of the shape of the diffuse X-ray emission from 
Wd 1. The diffuse X-ray emission shows an ellipticity of 0.75 
(corresponding to an eccentricity of 0.25) with a major axis 
along PA = 13° + 3° and a characteristic radius (half-width at 
half-maximum) of 25" (corresponding to 0.44 pc for a distance 
of 3.55 kpc). Muno et al. ( 12006b ) discuss possible origins of 
the diffuse X-ray component. While the emission in the core of 
Wd 1 can be attributed to colliding winds of O and W-R stars 
with initial masses >15 Mq, less than 30% of the extended dif- 
fuse emission can be explained by a spatially dispersed pre- 
main sequence star population with masses <2 Mq. A possible 
explanation for a;70% of the extended component of the dif- 
fuse emission are nonthermal particles accelerated by colliding 
winds and supernova remnants. 

The fact that both the amount of ellipticity and the direction 
of the major axis of the diffuse X-ray emission are in good 
agreement with the structural parameters derived from the NIR 
data for stars with masses between 10 and 30 M© supports the 
idea that colliding winds from the most massive stars give rise 
to the diffuse X-ray emission in the core of Wd 1 . 

In their study of the structure and dynamics of the ONC, 
Hillenbrand & Hartmann (1998) find an eccentricity of 0.30 for 
ONC. They conclude that rotational flattening of ONC seems 
unlikely and that the elongation might reflect the shape of the 
parental molecular cloud. 

Given the still young age of Wd 1 of 3 to 5 Myr, its elon- 
gation, too, should trace back to the formation of the clus- 
ter Possible explanations are formation out of an elongated 
parental molecular cloud, formation of the cluster out of two 
or multiple molecular sub-cores, out of a rotating molecular 
cloud, or a combination of these. 

Table 8. The decreasing half-mass radius with increasing stel- 
lar mass gives evidence for the presence of mass segregation in 
Wd 1. The value in the mass bin 30^0 Mq is from C05. 



mass bin 

fhrn 



3-10 Me 
1.1 pc 



10-30 Mo 



30-40 Mo 



0.8 pc 0.44 pc (25") 



5.2. Cluster dynamics 



According to Binney & Tremaine (I1987I I. the velocity disper- 
sion cr of a relaxed stellar system in virial equilibrium is 

0-= VO.4 X G X Mci/rhm 

For a cluster mass of > 45,000 M© and a 2-D half-mass ra- 
dius of rhm = l.Opc (see Fig. [8] right), we derive a 1-D veloc- 
ity dispersion cr > 4.5km/s. This is in good agreement with 
the velocity dispersion of 5.0 + 1.7km/s derived by Mengel & 
Tacconi-Garman (2007) based on a weighted average for ra- 
dial velocity measurements of four red supergiants in Wd 1 . In 
terms of proper motion, 4.5km/s corresponds to 0.25 mas/yr 
at a distance of 3.55 kpc. As shown by Stolte et al. (120071) for 
the Arches cluster, multi-epoch adaptive optics observations at 



8 to 10 m class telescopes now achieve measurement accura- 
cies better than 0.5 mas/yr. Thus it should be possible to mea- 
sure the velocity dispersion of a significant number of cluster 
stars using two epochs of high-angular resolution observations 
separated by a few years. From this, any deviation from virial 
equilibrium due to the cluster's young age can be investigated. 
Rotation of the cluster as a whole might result in an anisotropy 
in the observed velocity dispersion of stars. While in the di- 
rection parallel to the rotation axis stars should just exhibit the 
internal velocity dispersion of the cluster, in the direction per- 
pendicular to the rotation axis stars should on average show a 
larger velocity. 

The crossing time for a star located at the half-mass radius 
and with a velocity equal to the velocity dispersion is tcross = 
I'hm/cr = 3x10^ yr. For a cluster age of 3.6 Myr, this corresponds 
to a! 12 crossings since the formation of the cluster, indicating 
that the cluster already has experienced significant dynamical 
relaxation and mass segregation. 

Assuming a total number of cluster members Ntot - 10^ 
(Table |6]), we can also compute the half-mass relaxation time 
(see, e.g., Spitzer & Hart 119711 ): 



0.2N„ 



^relax ~ 



ln(0.1N,o,) 



Va4xt„oss=4xl0*'yr 



Thus even if the gaseous and stellar mass of the cluster 
would have been constant since its formation, it would not yet 
be dynamically relaxed, and further mass segregation would be 
expected. The expulsion of any remnant gas tracing back to the 
formation of the cluster by the multiple supernovae explosions, 
as well as the loss of 12% («i6,000Mo) of the total initial steflar 
mass of the cluster since its formation due to stellar winds and 
supernova explosions keep Wd 1 out of virial equilibrium. 

The effect of rotation on the dynamical evolution of a clus- 
ter was discussed by Ernst et al. (2007). They found that the 
presence of rotation accelerates the dynamical evolution only 
in systems with a single mass stellar component. In clusters 
with a realistic mass function, this effect disappears as mass 
segregation becomes the primary means to accelerate core col- 
lapse of the cluster Thus we do not expect that any rotation of 
Wd 1 should have a strong effect on its dynamical evolution. 

As discussed, e.g., by Portegies Zwart et al. ( 12007b . the 
presence and fraction of binary stars in a cluster can also have a 
strong effect on its dynamical evolution. Crowther et al. (I2006I I 
discuss evidence that at least between 1/3 to 2/3 of the W-R 
stars in Wd 1 are members of binary systems. Bonanos (I2007I I 
identifies 5 eclipsing binaries among the massive stars in Wd 1 . 
A systematic survey for binaries and multiple systems in Wd 1, 
in particular among the intermediate to low mass stars, is still 
missing. Thus we currently cannot quantify the binary fraction, 
and its effect on the dynamical evolution of Wd 1 . 

6. Comparison to the NGC3603, Arches and 30 
Doradus starburst clusters 

In Table l9] we present a compilation and comparison of astro- 
physical quantities of the galactic young massive clusters Wd 1 
with ONC, NGC 3603 YC, and Arches, as wefl as the R136 
cluster in the 30 Doradus region in the Large Magellanic Cloud. 
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Table 9. Comparison of Westerlund I's astrophysical parameters with other young, massive clusters 





Wdl 


ONC 


Wd2'' 


NGC 3603Ye' 


Arches 


R136 cluster'' 


Age [Myr] 


3.6 ±0.7 


0.3-1.0 


1-3 


1-3 


2-3 


2-3 


trekx [Myr] 


400 


2.1 






2? 




Tcore [pc] 


0.4 


0.16-0.21 


0.002-0.2 


7 


0.034-0.30 




fhm [pc] 


1.0 


0.8 




0.002-0.4 


0.24 


1.1 


e 


0.15-0.20 


0.30 


? 


? 


? 


? 


N,„, 


10' 


>2800 




> 1680 


>3450 




M,„, [Mo] 


5x 10" 


0.18 X 10* 


0.7 X 10* 


> 0.7 X 10* 


> 1.3 X 10* 


10 X 10* 


(mass range in Mq) 


0.4-120 


0.1-35 


0.08-120 


0.1-120 


1-120 


0.5-25 



" Hillenbrand & Hartmann ( 1998 1, and references therein. 

*" Ascenso et al. (20071, and references therein. 

■^ Stolte et al. ( 2005 2 006b , and references therein. 

■^ Figer et al. (1999), Stolte et al. ('2005^, and references therein. 

° Andersen et al. C2007j , Brandl et al. (.1996 J . and references therein. 



All clusters have in common that they posses at least several 
O-type stars with masses in excess of 30 Mq. The large num- 
ber of evolved massive stars, including hypergiants and W-R 
stars, present in Wd 1 gives evidence that it is the oldest, most 
evolved cluster in the sample. IVlass segregation seems to be 
a common feature among young clusters. As shown in Fig. l9] 
Wd I's half-mass radius increases with decreasing stellar mass. 
A similar trend has been observed by Brandl et al. d 19961 1 for 
the core radius of the R136 cluster Compared to NGC 3603 YC 
and the Arches cluster, Wd 1 has a rather large half-mass radius 
of ail. Ope, quite similar to the half-mass radii of the ONC and 
the R 136 cluster. 

Since dynamical mass estimates have not yet been derived 
for any of the clusters, and since the crowding and the presence 
of bright stars limits the ability to detect the low-mass stellar 
content of these clusters, mass estimates are in general lower 
limits. Wd 1 with a total stellar mass between 21,000 M© and 
52,000 Mq is the most massive, and with a total stellar popula- 
tion of up to 100,000 stars also the most populous among the 
galactic young massive clusters presented in Table l9] Only the 
R136 cluster in the 30 Doradus region in the Large Magellanic 
Cloud, for which Andersen et al. ( I2007I I derive a mass of 2.5 
to 3 xIO'^Mq counting stars with masses down to 2.35 M©, is 
more massive. 



7. Summary and Outlook 

We have analysed near-infrared NTT/SofI observations of the 
starburst cluster Westerlund 1, which is among the most mas- 
sive young clusters in the Milky Way. A comparison of colour- 
magnitude diagrams with theoretical main-sequence and pre- 
main sequence evolutionary tracks yields improved extinc- 
tion, distance and age estimates of Aks - 1.1 3 ±0.03 mag, d = 
3.55±0.17kpc (DM = 12.75±0.10mag) and t = 3.6±0.7Myr, 
respectively. We derive the slope of the stellar mass function for 
stars with masses between 3.4 and 27 Mq. In an annulus with 
radii between 0.75 and 1.1 pc from the cluster centre, we get a 
slope of r - -1.3, i.e. the Salpeter slope. Closer in, the mass 
function of Westerlund 1 is shallower with F = -0.6, while at 
larger separations from the cluster it is getting steeper, reach- 



ing F = -1.6 for separations larger than 1.5 pc. This is in good 
agreement with the change in mass function slope found in the 
starburst cluster NGC 3603 YC (Stolte et al.' 205g] l. 

Only considering stars with masses between 3.0 and 32 Mq, 
we derive a half-mass radius of 1 .0 pc for the cluster. Wd 1 ex- 
hibits clear deviations from spherical symmetry. The distribu- 
tion of stars with masses between 10 and 32 Mq has an eccen- 
tricity of 0.20 with the major axis aligned roughly in the North- 
South direction. The distribution of stars with masses between 
3.0 and 10 Mq is elongated in the same direction with an eccen- 
tricity of 0. 15. The flattening of the cluster might be explained 
by rotation of the cluster along its minor axis. 

By extrapolation of the observed mass function slopes for 
diff'erent annuli, we derive an upper limit of the total initial stel- 
lar cluster mass of a 52, OOOMq. Stellar evolution at the high 
mass end should have reduced the total initial cluster mass by 
~6,000 Mq over the past 4 Myr By adding up the individual ini- 
tial masses of stars directly detected in the cluster, we derived a 
lower limit for the total present cluster mass of m > 2 x 10"* Mq. 

With a present-day mass of 20,000 to 45,000 Mq and an 
initial stellar mass of a!52,000MQ, Wd 1 is the most massive 
starburst cluster identified in the Milky Way to date, and about 
10 times as massive as the ONC, and 2 to 4 times as massive a 
NGC 3603 YC. Additional cluster members located outside the 
field of view of the SofI data might push the mass of Wd 1 even 
higher 

In a following paper, we will report of the analysis of high- 
angular resolution adaptive optics data, which trace the pre- 
main sequence population of the cluster down to aO.2 Mq. 

Further studies should aim at determining proper motions 
for individual cluster members as well as more radial velocity 
measurements. From this precise models of the cluster kine- 
matics and dynamics could be derived, which also might enable 
one to trace back the dynamical evolution of the cluster since 
its formation a!4Myr ago. Wd 1, being the most massive, and 
one of the most nearby starburst clusters, could turn into the 
prime example for studies on starburst cluster formation and 
evolution. 
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